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ABSTRACT 

We present photometric and spectroscopic results for the low mass pre-main 
sequence (PMS) stars with spectral types K - M in the cometary globule (CG) 
30/31/38 complex. We obtained multi-object high resolution spectra for the 
targets selected as possible PMS stars from mult i- wavelength photometry. We 
identified 11 PMS stars brighter than V = 16.5 with ages < 5 Myr at a distance 
of approximately 200 pc. The spatial distribution of the PMS stars, CG clouds, 
and ionizing sources (O stars and supernova remnants) suggests a possible trig- 
gered origin of the star formation in this region. We confirm the youth of the 
photometrically selected PMS stars using the lithium abundances. The radial 
velocities of the low mass PMS stars are consistent with those of the cometary 
globules. Most of the PMS stars show weak Ha emission with Wx{Ha) < 10 A. 
Only 1 out of the 11 PMS stars shows a moderate near-IR excess, which suggests 
a short survival time (t < 5 Myr) of circumstellar disks in this star forming en- 
vironment. In addition, we find five young late type stars and one Ae star which 
have no obvious relation to the CG 30/31/38 complex. We also discuss a possible 
scenario of star formation history in the CG 30/31/38 region. 
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1. Introduction 

Most massive stars are found in their natal clusters and associations. Low mass pre- 
main sequence (PMS) stars are also found in close proximity to massive stars, such as seen 
in Orion Nebula Cluster (O'Dell et al. 1993; Bally et al. 1998), and in vicinity of the Wolf- 

Rayet + O star binary system, 7^ Vel (Pozzo et al. 2000). As soon as the massive stars begin 
to burn hydrogen, the protostellar clouds and T Tauri stars in the vicinity of the massive 
stars could be affected by the winds and strong radiation of massive stars. In contrast to 
stars formed in such a strong radiation field environment, low mass stars formed in a quiet 
environment unaffected by massive stars, such as in the Taurus star forming region, may be 
able to retain larger reservoirs of circumstellar materials, which could prolong the lifetime 
of accretion disks. 

Globules are isolated dense clouds of molecular gas (Bok & Reilly 1947) with masses of 
order of 10 - 100 Mq and typical sizes of 0.1 - 1 pc. Bok & Reilly (1947) and Bok (1977) 
initiated studies of star formation in globules, and these have been followed by many others 
(Schwartz 1977; Reipurth 1983; Sugitani et al. 1986, 1995; Lefloch & Lazaref 1994; Megeath 
et al. 1996; Hester et al. 1996; Alves et al. 1997; Wolk & Walter 1998; Bally et al. 2001; Kim 
2002; Kim et al. 2003, and references therein). These globules may form one to a few low 
mass stars in a very short time scale compared to the mean life time of globules, 1 - few x 10^ 
years, (see Elmergreen 1998 for review; and references therein). These globules are possible 
birth places for isolated low mass stars, such as the Sun. 

In this paper, we present a study of the star forming region CG 30/31/38 , which lies 
about 50 - 100 pc from an early O star and OB associations. This distance range is still close 
enough to be influenced by the radiation, but too far to be affected by the winds of the hot 
stars. The CG 30/31/38 star forming region is a complex of small molecular cloud cores in 
the Gum Nebula evaporating due to ionizing radiation from an O star, ( Pup, and perhaps 
also due to progenitor of the Vela supernova remnant (SNR). The total ionizing radiation of 
two OB associations. Vela 0B2 and Trumpler 10, may also be responsible for some ionizing 
fluxes in the Gum Nebula, and might have affected in this CG 30/31/38 region. 
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1.1. The Gum Nebula 

The Gum Nebula, one of the largest H II regions in our Galaxy, was discovered about 
half a century ago (Gum 1952, 1955). It covers a large part of southern sky, in Vela, Puppis, 
Pyxis, Canis Major, and Carina. It is located in the Galactic plane, with a 36° diameter, 
centered at I — 258° and b — —5°. The nebula is the largest apparent H II region, with a 
36° diameter. The hnear diameter is 250 pc at a distance of 450 pc (Brandt et al. 1971; 
Reynolds 1976a,b). The nature of the Gum Nebula remains unclear: this may be an old 
supernova remnant (SNR), a fossil Stromgren sphere (Brandt et al. 1971), expanding H II 
regions (Reynolds 1976a,b; Chanot & Sivan 1983) and/or a stellar wind bubble (Sahu & 
Sahu 1993, and references therein). The hot stars ( Pup (04If) and 7^ Vel (WG8+07.5I), 
and two OB associations (Tr 10 and Vela 0B2) near its geometric center can account for 
most of the photoionization in the nebula, and their UV radiation probably is the cause of 
the evaporation of the CGs (Bertoldi & McKee 1990; Lefloch & Lazaref 1994). 

The distances toward these associations, ( Pup, 7^ Vel, and the Vela SNR, and their re- 
lation with the Gum Nebula are poorly understood, and remain controversial. The canonical 
values for the distances to these objects are about 400 - 500 pc (Brandt et al. 1971; Hawar- 
den & Brand 1976; Zealey et al. 1983; Henning & Launhardt 1998, and references therein). 
However, more recent studies, especially using the Hipparcos data suggest distances of 200 - 
400 pc (de Zeeuw et al. 1999; Knude & Nielsen 2000; Hoogerwerf et al. 2001, and references 
therein) . 

The Gum Nebula contains at least 32 cometary globules (CGs) (Hawarden & Brand 
1976; Sandqvist 1976; Reipurth 1983; Zealey et al. 1983). Cometary globules are small bright 
rimmed isolated clouds with head and tail cometary morphology. These are evaporating 
clouds which are normally associated with H II regions and OB stars. High resolution 
imaging with the Hubble Space Telescope has revealed smaller, solar-system size globules in 
Orion (O'Dell et al. 1993; Bally et al. 1998) in which dusty protostellar disks are seen in 
silhouette. Herbig-Haro (HH) objects and outflow sources are often found in Bok globules, 
cometary globules (CGs), and bright rimmed clouds. Examples include B 335 (Frerking & 
Langer 1982; Keene et al. 1983), HH 120 in CG 30, HH 46/47 in the Gum Dark Cloud (GDC) 
1-7 (Bok 1977; Reipurth 1983), and CB 34 (Keene et al. 1983; Clemens & Barvainis 1988; 
Khanzadyan et al. 2002). Some PMS stars have been identified in vicinity of Bok globules 
and CGs (Reipurth 1983; Pettersson & Reipurth 1994; Alves et al. 1997; Kim 2002). 

The tails of the CGs in the Gum Nebula are directed away from the central ionizing 
sources, ( Pup, 7^ Vel, and Vela pulsar, suggesting that CGs are evaporating away from 
a common center. The locations of CGs and the direction of the tails show that CGs 
probably he at the periphery of the wind-blown shell. If CGs were located on the hne of 
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sight between us and their ionizing sources, their tails would not be noticeable. However, the 
CGs located along the periphery of a sphere can be seen easily, because the tails evaporate 
in the transverse direction as projected in the sky. Therefore, even though CGs are located 
at about 50 - 70 pc away from the ionizing sources, the dispersion in distances among the 
CGs can not be too large, for we are most likely to sample those near or along the periphery 
that show prominent tails. 

From his kinematic studies of the CGs in the Gum Nebula, Sridharan (1992) showed 
that CGs are expanding away from the common center with a velocity of ~12 km s~^. The 
velocity gradients along the CG tails imply that these tails have expanded for about 3 - 
6 Myrs. An extended shell around Vela 0B2 association, the IRAS- Vela Shell, has a size 
and location that seems to be consistent with the ring-like spatial distribution of the CGs 
(Sahu 1992). Note that the radius of this ring of CGs are smaller than that of the Gum 
Nebula. The expansion of this shell might have been caused by radiation driven stellar winds 
and multiple supernova explosions. The expansion velocities of CGs and this shell, on the 
order of 10 km s~^, are consistent with this, implying that the cause of CG expansion may 
be the same as the one for the IRAS- Vela shell. 



1.2. The CG 30/31/38 complex 

Cometary globules, such as the CG 30/31/38 complex in the Gum Nebula, provide a 
unique and convenient setting wherein to study star formation in an intermediate radiation 
field at a moderate distance from hot stars. Figure 1 is the Digital Sky Survey (DSS) image 
of the CG 30/31/38 complex. CG 38 is the smallest CG in this complex, seen just below 
CG 30. The head of CG 30 (~2' in diameter) contains a Herbig-Haro system, HH 120, and 
the IR source CG 30-IRS 4 (Reipurth 1983; Pettersson 1984). Its ~10' long tail points away 
from the central ionizing sources (i.e. Vela SNR, Vela 0B2). CG 31 has 5 distinct clouds. 
The tail of CG31 A is ~25' - 30' long (1 - 1.5 pc at 200 pc). The change in the orientation 
of the tails of CG 31 coincides with the proper motion direction of ( Pup. Pettersson (1987) 
identified three early M stars with Wa(Hq,)>10 A in this region as Ha emitting T Tauri 
stars. 

Recent studies of the reddening in this direction (Nielsen et al. 2000; Knude & Nielsen 
2000) suggest a distance as close as 200 pc, about half the canonical value. We use 200 pc 
as a conservative value for the distance of this cometary globule complex throughout this 
paper. 
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2. Observations 

To identify PMS stars in the CG 30/31/38 region, we obtained X-ray images, opti- 
cal, and near IR photometry. We then obtained high resolution multi-object spectra of 
photometrically-selected samples to confirm the youth and membership of the samples. In 
this section, we describe these observations, the data reduction processes, and the calibration 
of the multi-wavelength photometry and high resolution multi-object spectroscopy. 

2.1. ROSAT/HRI photometry 

X-rays are an efficient means of selecting candidate low mass PMS stars (Montmerle 
et al. 1983; Walter et al. 1988, and references therein), because these stars are rapidly 
rotating and highly convective, the two ingredients which generate strong solar-like-coronal 
and chromospheric activity. Magnetically active PMS stars (G, K, and M dwarfs) have 
X-ray luminosities of >10^^ erg/sec, and X-ray to bolometric flux ratio {fx/ hoi) of about 
10-3 - 10"^ (Walter et al. 1988; Feigelson et al. 2002; Flaccomio et al. 2003a, and references 
therein). In contrast the Sun's X-ray luminosity varies from 10^^ to 10^ erg/sec, and its 
flux ratio varies from 10^^ - 10"'^. We use the spatial coincidence of bright stars with an 
X-ray-to-bolometric flux ratio (fx/fboi) > 10"^ as the initial selection of the potential PMS 
stars. Since this criterion also selects magnetically active foreground stars and fast rotators, 
we also obtained optical and near IR imaging data to select an unbiased photometric sample 
of candidate PMS stars. 

The CG 30/31/38 region was targeted in the ROSAT observation RH 201400, a 18364 
second exposure obtained between 12 and 15 May 1994. The image was centered at a=8^^9.2"\ 
5=-36°6 . Thirteen X -ray sources were detected in the 23' radius field of CG 30/31/38 re- 
gion by the standard SASS source detection algorithm in 12"xl2" or 24"x24" detection 
cells (Table 1). Eight of the nine sources detected in the 12"xl2" were recovered in the 
larger detection cell. These X-ray sources show a non-uniform spatial distribution, generally 
outlining the heads of the CG 30 and CG 31 (Figure 1). 

The ROSAT high-resolution imager (HRI) has little intrinsic energy resolution; we con- 
verted the HRI count rates (0.4 - 2 keV) to fluxes using the standard energy-to-counts 
conversion factor (ECF)^ for Raymond-Smith spectral model. We assume a 1 kT thermal 
spectrum, and converted Ay to Nh (Nh=2.0x10^^ cm~^ mag~^ x Ay). We estimated Ay 



^scc Table 8c and Figure 30 from The ROSAT High Resolution Imager (HRI) Calibration Report, 
ftp: / / sao-ftp.harvard.edu/pub/rosat/docs/hri_report99/hri_99.ps 



-6- 



values from the spectral types (see section 3.3) and observed B — V colors, adopting intrinsic 
colors from dwarf stars (Kenyon & Harmann 1995; Luhman 1999; Leggett et al. 2001). We 
assumed the standard galactic extinction law with R=3.1. 

2.2. Optical followup photometry 

We obtained optical photometry of selected CG complexes in the Gum Nebula on the 
nights of January 29 - 30, 1996. We used the CTIO 0.9m telescope with the 2048x2048 
Tektronix CCD (a pixel scale of O'.'4/pixel in a 13!6 field of view) and UBVRI filters. During 
March 6 - 12, 2001 we obtained more BVRI photometry in these fields with the same setup. 
All data were obtained on photometric nights, with the mean seeing of 1'.'2. 

Bias, dome flat flelds, and sky flat flelds were taken at the beginning and the end of 
each night. Long (60 and 300 sec) and short (10 and 20 sec) exposures were taken for all the 
object fields to avoid saturation of bright stars in the fields. Landolt (1992) standard fields 
were observed several times per night to establish photometric calibration. 

Since the Gum Nebula lies in the Galactic plane there are typically more than 2000 
stars in each field, many of which are heavily blended. We reduced the images with 
IRAF/DAOPHOT package. Since the point spread function (PSF) depends on the loca- 
tion within the CCD, we selected between 40 and 80 PSF stars in each field, and determined 
the spatial dependence of the PSF using quadratic fits in both directions of the CCD. 

The full width half maximum values used for the IRAF DAOPHOT photometry were 
between 1'.'2 to I'.'S pixels, while the radius for standard star photometry was 18 pixels. 
Therefore, an aperture correction was made to the 18 pixel (7'.'2) aperture. The photometric 
errors are about 1% at ^ < 18, and 5 - 10% at V ~ 21. The aperture correction is accurate 
to 1 - 2% in general. 

We used imwcs in wcstools (Mink 1999) version 2.7.2 for the astrometric solutions. 
The positions of stars were determined using the Gaussian centroid algorithm in IRAF phot 
routine, then the astrometric solution was derived by fitting the pixel coordinates of stars 
to the known positions of USNO-A2.0 astrometric standards (Monet 1998) located in each 
image. The imwcs package projects the image onto the tangent plane of the sky, rotates the 
image, and fits a polynomial to the x and y coordinates. In each image, about 120 - 150 
stars were matched with USNO-A2.0 astrometric standards with an average error of ~0'.'l. 
We then averaged the positions measured in the V, R, and I band images for each star to 
give the final position. 
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2.3. Neeir-IR photometry 

We obtained JHK images covering all the X-ray sources using the CIRIM detector on 
the CTIO 1.5m telescope during February 5 - 6, 1996 and February 17 - 20, 1997. CIRIM 
uses a 256x256 HgCdTe NICMOS 3 array, and at f/13.5 gives a pixel scale of O'.'44/pixel. 
Images were dithered by ±15" in RA and DEC, and median filtered to determine level of 
the sky. For the targets we used a 6-point dither pattern. We used a 4-point dither pattern 
for the standard stars. 

We reduced and analyzed the data using the IDL/DoCIRIM sottwaie (Walter 2000) and 
custom written IDL codes. We first generated median sky filtered and subtracted images. 
Standard stars were chosen from the Ehas (Ehas 1982) and UKIRT FS standard star hsts.^ 
For astrometry we match stars in an image with USNO A2.0 catalog stars, and determine 

the image center and plate scale. We performed standard star photometry using a radius 
of 20 pixels (8'.'8), and generated the photometric solution on the CIT system. For object 
stars we used a default aperture radius of S'.'l, but we also used smaller apertures for objects 
that are faint or close to neighboring stars. We derived the aperture corrections from the 
standard stars, and corrected our targets to 20 pixel apertures. In this study, we used only 
those stars with photometric error <5%, which corresponds to K ~15 - 16 magnitudes. 

2.4. Matching optical, neeir-IR counterpEirts, and X-ray sources 

We matched the coordinates of X-ray sources with those of optical and near-IR stars 
to identify potential counterparts. The ROSAT/HRI X-ray source positions have data has 
much larger positional error (see Table 1) than those of the optical and near-IR catalogs. 
The matching algorithm, therefore, identified optical and near-IR sources within the X- 
ray error circles (Table 1). The HRI flux hmit was such that a star matching the fx/ fboi 
>10~^ criterion has V < 16. Since X-ray, optical, and near-IR data were not obtained 
simultaneously, and PMS stars are variable, we also considered stars up to 1 mag fainter. 
Twelve of the HRI sources have bright (V < 17; K < 13) optical and near-IR counterparts 
within the 3a X-ray error circles. For the X-ray sources that were matched with more than 
one optical and IR source, we chose the brightest sources among all the sources within the 
error circle. The faintest optical counterpart, XRS 3, has V magnitude of 17.27; XRS 1 does 
not have an optical counterpart with V < 17 within its 3a X-ray error circle. 



^http://www.jach.hawaii.edu/UKIRT/astrononiy/calib/photjcal/fs -fundamental. html 
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2.5. Optical spectroscopy 

During March 6 - 8, 2002 we obtained multi-object high resolution spectra of the can- 
didate PMS stars (see Section 3), using the HYDRA multi- fiber spectrograph on the CTIO 
Blanco 4m telescope in echelle mode. The goals of these spectroscopic observations were to 
confirm the youth of these stars using Li I 6708A absorption lines and to measure their radial 
velocities. The spectra also reveal the spectral types, and we can estimate extinction toward 
the PMS stars by comparing the measured colors with the intrinsic colors expected for the 
observed spectral types. The radial velocities let us determine membership probabilities, by 
comparing the distribution of velocities with that of the CGs. Using CTIO/HYDRA, one 
can obtain spectra of up to 135 objects (stars and sky) in a 40' field of view simultaneously. 
Since the overall efficiency of the HYDRA in echelle mode is about 3%, we chose targets that 
were brighter than V ~ 17 mag. The S/N was sufficient to detect Li I lines with strengths 
expected of PMS stars. 

2.5.1. Spectroscopic sample selection 

We selected as targets those stars that are coincident with the HRI X-ray sources and 
other stars whose colors lie in the appropriate region in the Color-Magnitude Diagram (CMD) 
for PMS stars at a distance of ~200 pc. V magnitudes of targets range between 10 and 
17. In each field there are between 50 and 200 objects near the PMS locus in the CMD 
(Figure 2). At the bluer end {V — I ~1.0) of the CMD, the non- X-ray emitting stars in the 
PMS locus are heavily contaminated by background stars in the galaxy, but the fainter and 
redder candidates are less likely to be contaminated. Since the Gum Nebula is located in 
the Galactic plane, some of the fields suffer severe contamination. 

2.5.2. Spectroscopic observations, data reduction, and wavelength calibration 

We observed 144 objects in two 40' diameter fields on March 6 and 7, 2002. Three 
exposures were taken of each target field. These were combined using the IRAF combine 
with the median sigma clipping algorithm. Three offset sky spectra were also taken at a 
position offset 10" from the targets after each target exposure. These sky spectra have 10% 
of the integration time of the target spectra. These sky offset images were used to identify 
background nebula emission line profiles. The sky subtraction is discussed in more detail in 
section 2.5.3. 

We used the IRAF/HYDRA package for basic data reduction, and also some custom 
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written IDL codes for various calibration purposes. We first fit the CCD overscan of the each 
image interactively, and trimmed all images. We then created a combined bias (zero) for 
each night, and subtracted it from all images. We created a milky flat image from daylight 
sky exposures, which we used to flatten all other images. 

Scattered hght was subtracted from images using abscat task. The scattered hght is a 
global 2-D feature, and is subtracted before each aperture is traced and extracted. Then the 
IRAF/apall task was used to extract spectra. The FWHM of a fiber is about 6-7 pixels, 
so we selected 3 pixels for upper and lower widths, and extracted a total of 6 pixels centered 
at the center of each trace. 

We took about 5 dusk and 5 dawn twilight sky spectra each night to make combined 
dusk and dawn sky spectra. Th-Ar comparison lamp exposures were taken at the beginning 
and end of each night for wavelength calibration. The mean residual of the wavelength 
solutions is < O.OSA (0.5 pixel). Etalon comparison lamp exposures were taken before and 
after taking Th-Ar lamps, and between every configuration of object fields to correct for the 
zero point shift throughout the night. We selected as reference the etalon spectrum which 
was taken immediately before or after the Th-Ar comparison exposure each night; all other 
etalon spectra were cross-correlated in pixel space. We calculated the zero point shifts in 
pixels, and then applied the shift to derive wavelength vectors for all object spectra. 

Since the Th-Ar signal was weak during the first night of observation, because of a 
faulty lamp, we used the second night's Th-Ar solution. Some apertures that contained 
high enough S/N in first night's Th-Ar spectra were cross-correlated with the second night's 
Th-Ar spectra. The median pixel shift of 0.023 pixels is very small. We applied this shift to 

first night's Th-Ar solution; the rest of the wavelength calibration using the etalon was done 
the same way as on night 2. IRAF/dispcor was then used to dispersion-correct sky spectra 
and object spectra using Th-Ar solution as reference spectrum. The starting wavelength and 
dispersion (A/pix) were derived for each of the 135 fibers; then the zero point was shifted 
using the cross-correlation results from Th-Ar and etalon spectra. 

2.5.3. Sky Background Subtraction 

The nebulosity in these regions complicates the sky subtraction. In addition to the 
offset sky spectra, we also devoted between 10 and 20 fibers in each configuration field to 
sky. Here we use the term sky background to refer to the emission lines due to nebulosity, 
which is particularly important for measuring the strength of the stellar Ha lines. The sky 
emission fine strengths of Ha and [N H] vary across the 40 ' field of view. The background 
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sky line profiles vary spatially due to the kinematics of the expanding gas shell in the Gum 
Nebula. 

We group the sky spectra into three general categories, based on the Ha and [N II] hne 
profiles and strength. Group 1 has narrow and strong Ha and [N II]; group 2 has asymmetric 
and broader line widths than group 1; and group 3 has broader line widths than group 2, 
possibly with double-peaked lines. 

To check the strength and profile of night sky emission lines, we compared the night 
sky spectra obtained during the long exposure with the combined three offset sky spectra 
which were taken 10" offset from the target field immediately following each of the three 
sub-exposures. Although the strength of the Ha lines in the offset sky spectra were 10% 
of long exposure sky lines, the offset sky spectra helped determining which subgroup of sky 
profile should be used for each star. 

We subtracted the appropriate sky spectrum from the object spectra. Due to the spa- 
tially variable background sky line profiles, sky background can not be perfectly subtracted. 
However after careful checking, we are confident that we can derive the equivalent width of 
Ha to an accuracy of ~0.8A. Line emission is not a problem at the A6708 A wavelength of 
Li I. 



3. Results 

From X-ray, optical, ncar-IR photometry, and optical spectroscopy, we have identified 
16 young stars (Table 2). Eight of the 13 X-ray sources (XRS) listed in the Table 1 arc 
spectroscopically confirmed young stars, two X-ray sources were not observed, and three 
do not appear to be young stars. We identified eight more young stars, not detected in 
the ROSAT/HRI observations, based on the Li I A6708 A absorption line. Of the 16 stars, 
radial velocities of eleven are consistent with that of the CGs (Sridharan 1992; Nielsen et al. 
1998). In the next section we discuss our results from X-ray photometry; optical and near 
IR photometry; and optical spectroscopy. 

3.1. X-ray photometry 

In Table 1 we present X-ray count rates of the 13 X-ray sources. We also present the 
offsets in X-ray coordinates and the optical counterparts. Table 2 lists spectroscopically 
confirmed young stars along with their spectral types. Ha and Li equivalent widths, Li 
abundances, and Ay values. 
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Using the ROSAT/HRI photometry, optical photometry and spectroscopy, we derived 
X-ray to bolometric flux ratio for 10 X-ray sources (Table 1). Using the spectral types (see 
section 3.3.1), we adopted bolometric correction values from a color table that was compiled 
using published color tables by Kenyon & Harmann (1995); Luhman (1999); Leggett et al. 
(2001), which can be found in Table 3 of Sherry et al. (2004). 

We used Table 3 from Sherry et al. (2004) to estimate unreddened colors and bolometric 
corrections for our spectroscopic sample. This table combines spectral types, bolometric 
corrections, and colors from Kenyon & Harmann (1995) and Leggett et al. (2001) with the 
PMS temperature scale of Luhman (1999). We derived Ay for each star by comparing its 
predicted unrcddenned B - V color to the observed B - V color. As a check on the derived 
AyS, we compared the spectral energy distribution (SED) for each star derived from our 
V, R, I, J, and H band data to model SEDs for K and M stars with a range of AyS. AyS 
derived from the SEDs were consistent with the Ays we calculated from the B - V colors, 
but with larger uncertainties on Ay. This suggested that our uncertainties on Ay were two 
- three times as large as we had estimated from the B - V colors. The larger uncertainties 
are incorporated into Table 2 and Figures 3 and 4. 

3.1.1. fx/fv and fx/ hoi 

Figure 3 shows X-ray flux vs. bolometric flux ratio of 10 X-ray sources listed in Table 1 
in logarithmic scale (log(/x) vs. \og{fhoi))- XRS 5 and XRS 12 are confirmed not to be PMS 
stars, and since we do not have optical flux of XRS 9, we do not include these three sources 
in this plot. Two high fx points are KWW 1863 (XRS 4) and KWW 1637 (XRS 6). These 
high count rates suggest strong coronal activity or flaring. All the sources appear above 
logifx/fboi) = -3 line. 

Figure 4 is a X-ray to bolometric flux ratio vs. effective temperature in logarithmic 
scale {log{fx/ fboi) vs. log(Tejj) plot of the 10 X-ray sources plotted in Figure 3. The 
\og{fx/ fboi) values ranges from ~-3 up to ~-1.5. Active K and M type stars arc known to 
have log{fx / fboi) ratio from —4 up to —1, while log{fx/ fboi) of PMS stars with spectral 
types earlier than G type range between -4 to -2. 

In Figure 4 we overplot 3a detection limits. These stars lie just at or slightly above 
this 3(7 detection limit. The mean log{fx/fboi) for other star forming regions is about —3.5 
(Feigelson et al. 2002; Flaccomio et al. 2003a,b, and references therein). This suggests that 
the HRI sources represent <15% of the X-ray bright PMS population if we assume that 
Lx/Lboi of the 0.1 — 2 Mq PMS stars follows that of typical young star clusters. If the 
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log{fx/fboi) of this star forming region follows the distribution of ONC (Flaccomio et al. 
2003a,b), for the mass range of 0.2 Mq - 1.0 M©, there could be roughly six times more 
PMS stars with log{fx/fboi) < -3.0. 

3.1.2. Variability 

We examined the X-ray light curves for variability. Since the number of source counts 
is small, we could look only for evidence of gross changes in the count rate. Of the brighter 
sources, only KWW 1637 (XRS 6), a double-lined spectroscopic binary PMS system, shows 
any evidence of variability, at about the 2a confidence level, in 30 ksec bins. The variation 
on timcscalcs between 30 and 250 ksec is no more than about 30% from the mean. Among 
the weaker sources, KWW 1055 (XRS 10) appears variable. All the counts arrived in the 
first half of the observation. The probability that the source is constant is 0.007 based on 
statistics. 

3.2. Color-magnitude diagram 

Figure 2 is the V,V — I CMD for objects in the CG 30/31/38 region. Stars that are 
likely counterparts of X-ray sources are indicated with crosses and asterisks. The limiting 
magnitude of X-ray flux, assuming \og{fx/ hoi) = -4, corresponds to V magnitude of ~16. 
The reddening vector, plotted on upper part of the CMD, is a line running almost parallel to 
the isochrones. Evolutionary tracks and isochrones are from Bareffe et al. (1998) modified 
using temperature scale and colors of Kenyon & Harmann (1995); Luhman (1999); Leggett 
et al. (2001) for low mass stars. Isochrones for ages of 2 Myr, 5 Myr, and 2 Gyr are plotted 
for a distance of 200 pc. 

The positions X-ray sources plotted in the CMD (Figure 2) seem to delineate two distinct 
loci parallel to either the ZAMS or a PMS isochrone. The upper locus is defined by six X-ray 

sources (x symbols) and three previously cataloged PMS stars (diamond symbols; Herbig 
& Bell (1988)). A 5 Myr isochrone (Bareffe et al. 1998) at the 200 pc distance of the cloud 
complex is a good match to this locus (see Figure 2). Four X-ray sources (asterisks) and one 
Ha source (diamond; Schwartz et al. (1990)) trace a distinctly different locus about 2-3 
magnitudes fainter than the upper locus. At a 200 pc distance these objects lie below the 
main sequence; they must represent a distinct and more distant population. 

By using the minimum likely distance of 200 pc to CG 30/31/38 complex (Knude et 
al. 1999), we can place an upper limit on the age of the stars if they are indeed PMS stars 
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at a common distance. The uncertainty in the age does correlate with the uncertainty in 
the distance. However, for young low mass stars one may use Lithium abundance to further 
constrain the age (sec §3.4). For example, M stars with log N(Li) ~3 are likely to be younger 
than ~5 Myrs old (Walter et al. 1994; Martin et al. 1994, and references therein). Conversely, 
by assuming (or measuring) an age we can constrain the distance to the cloud complex, if 
spatially related to the stars. 

3.3. Spectra 

Mult i- wavelength photometry enables us to produce a large and unbiased sample of 
candidate PMS stars, but is subject to contamination. Because photometry alone can be 
ambiguous, spectroscopy is necessary to confirm the youth and membership of stars using 
Li I A6708 A lines and radial velocities. In the previous sections, we described selection of 
candidate PMS stars using X-ray, optical, and near IR photometry. To confirm the youth of 
our photometrically selected PMS stars, we obtained high resolution multi-object spectra of 
our candidates. In section 2.5 we described the multi- fiber spectroscopic observations using 
CTI0-4m/HYDRA, data reduction, and calibration. Here we present the spectra and results 
from our spectroscopic observations and analysis. 

3.3.1. Spectral type determination 

Spectra of PMS stars are shown in Figures 5a - 5e. We fiattened the continua of the 
spectra using a polynomial fit, and normalized to unity. We estimated the spectral types by 
comparing with spectra of 26 standard stars obtained at comparable dispersions. 

The AA 6400 - 6720 A wavelength range is not ideal for spectral type determination, 
because of the paucity of strong lines that can be used for spectral type determination. We 
estimated the spectral types mainly using Fe I and Ca I lines. The strength and structure 
of the A6495 A blend was used to distinguish spectral types for spectral range of F - K. TiO 
bands were used for the determination of late K and M spectral types. 

The results are summarized in Table 2, and also discussed in the Appendix. The un- 
certainties on the spectral type are about 0.5 subtypes for K and M stars, and are larger for 
stars with earlier spectral types. 
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3.3.2. Spectroscopically confirmed PMS stars 

In Figures 5a - 5e wc show the spectra of the confirmed and possible PMS stars. Two 
of the confirmed PMS stars (KWW f892 and KWW f953) were observed twice. Three of 
the X-ray sources (XRS 5, 11, and 12) in the photometrically selected sample do not appear 
to be PMS stars. Spectra of KWW 1043 (XRS 8) and XRS 13 were not obtained, because 
these were too faint. We identified the Li I 6708 A absorption hne in 15 stars. Ten of these 
also have Ha in emission. Only KWW 873 shows near-IR excess. One star on the lower 
locus, KWW 314, is an early type star with an Ha profile reminiscent of an Ae star. The 
spectral type is A3 (see the Appendix). Not surprisingly, no Li absorption is evident. 

Based primarily on the presence of the Li I 6708 A line (next section) , we confirm that 
eight X-ray sources. KWW numbers 464, 1892, 598, 1863, 1637, 873, 1055, and XRS 9, are 
PMS stars. Ha line profiles of these PMS stars are described in more detail in the Appendix. 
We also find 8 other young or active stars that were not detected (or not observed) in X-rays, 
KWW 314, 975, 1125, 1302, 1333, 1806, 1953, and 2205 (see Table 2). 

Three PMS stars from the Herbig-Bell catalog (Herbig & Bell 1988), HBC 553, 554, 
and 555, also lie along the upper locus in the V, V — I CMD (diamonds in Figure 6). We 
observed two of these stars, KWW 1892 (HBC 553) and KWW 975 (HBC 554). We did 
not observe KWW 1043 (XRS 8), the hkely counterpart of HBC 555. Table 2 includes the 
derived properties of these HBC sources. 

Short descriptions of each of these stars can be found in Appendix. Table 2 summarizes 
our results. In Figure 6 the stars marked with x's and asterisks are the X-ray sources as 
in Figure 2. We mark the confirmed PMS stars with filled circles (upper locus) and open 
triangles (lower locus). Note that the PMS stars in the upper locus are clustered spatially 
below the CG 31 complex (boxes with circles in Figure 1). 

Inclusion of the optically-identified young stars strengthens the case for the existence 
of 2 parallel loci in this CMD (Figure 6) The upper locus is idefined by the X-ray sources 
(x), HBC sources (diamonds), and spectroscopically- young stars (filled circles). The Li 
abundance, fx/ hoi ratio, and the Ha emission, are all consistent with the interpretation of 
the upper locus as an isochrone defining a PMS stellar population. 

The lower locus is defined by X-ray sources (asterisks) and spectroscopically confirmed 
young stars (open triangles). Note that only one X-ray source (KWW 1055) is confirmed to 
be a PMS star, while XRS 11 and 12 are not PMS stars. As mentioned earlier, the spectrum 
of XRS 13 was not obtained. If the stars on the lower locus are at 400 pc, the age inferred 
from isochrone fitting is >100 Myr. If young, these stars are likely to be at a distance greater 
than 400 pc. They may be unrelated to the CG complex and the stars on the upper locus. 
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3.4. Lithium abundance 

The lithium abundance is of interest because it both probes stellar structure and age, 
and also constrains primordial nucleosynthesis models (Duncan 1991, and references therein). 
The universal primordial Li abundance, logN(Li), is 3.1 (on the scale of logN(H)=12), how- 
ever the Li abundance of halo stars is 2.1, while the mean distribution of logN(Li) values for 
PMS stars in young clusters peaks at 3.1 (Duncan 1991; Martin et al. 1994). 

High Li abundances can be indicative of the youth of solar mass or less massive stars, 
because the temperature at the base of convection zone is still too low to burn hthium for 1 
- 10 Myrs. These cool stars (M < 1 Mq) have deeper convective envelopes, and processes 
such as convective overshooting mix the surface and deeper (hotter) layers causing faster 
Li depletion than in higher mass stars. For higher mass stars (M > 1 Mq) Li depletion is 
not significant, because they spend less time in the PMS stage, and the convection zone is 
thinner than in low mass stars (King 1993; Ventura & Zappieri 1998). Because of this fast 
lithium depletion, lithium can be a good indicator of ages in low mass stars. 

We derive hthium abundances, log N(Li) using both NLTE (Pavlenko & Magazzu 1996) 
and LTE (Duncan 1991) analyses. The results are presented in Table 2. Since these curves 
of growth are valid for Teff > 3500 K, the curve of growth for lower temperature were 
extrapolated for the stars 3200 K < Teff < 3500 K. The derived Li abundances for K6 and 
K7 stars (XRS 6 and XRS 7) are higher than the universal value (3.82 and 3.56 for NLTE). 
The high value for these K stars are consistent with observations of other PMS stars, e.g., 
Basri et al. (1991); Martin et al. (1994); Walter et al. (1997). 

The individual values of log N(LI) among the M stars, all on the upper locus, range 
from 2.19 to 3.08 (in NLTE), with a mean of 2.68 ±0.28 (3.26 ±0.15 in LTE), where the un- 
certainties are the standard deviation of the means. The uncertainly on any single log N(Li) 
value is about 0.2 dex. The mean Li abundance is consistent with the expected age of the 
cloud complex (<3Myr). The wide range in log N(Li) among the individual stars could be 
due to the difficulty of determining the continuum level. The uncertainty in the adopted 
Teff of up to 200 K compounds the uncertainty. For KWW 1892 and KWW 975, the WA(Li) 
for two difi^erent nights vary by 0.14 and 0.16A, which caused 0.38 and 0.52 dex differences 
in the Li abundances. 

The F - G stars on the lower locus also show high Li abundances. Since F - G stars do 
not burn lithium as fast as late K or M stars do, such stars need not be pre-main sequence, 
but could be 50 - 100 Myr old main sequence stars. 
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3.5. Kinematics of PMS steirs and CGs 

Radial velocities were derived using Li I A6707.8 A and Ca I A6717.685 A lines. The 
FWHM of a fiber of the HYDRA/echelle spectrograph is about 6 pixels (0.96A, ~43.9 km s^^), 
however we can derive radial velocity good to about 3 km s~^ by centroiding. To estimate 
radial velocities, we derived 5X {Xobs ~ ^o), where Xobs is observed wavelength, and Aq for 
Li I and Ca I are from Rao & Lambert (1993). The heliocentric correction of ~-8.2 km 
was then applied. To compare the velocity with Ylsr (velocities with respect to the Local 
Standard of Rest) for the CGs, we then derived Ylsr for all the stars. The correction to 
Ylsr is about -17.3 km s~^ (Lang 1980). We present the radial velocities in Table 2. 

The Ylsr of the K - M type PMS stars range from 4.2 - 8.9 km s~^, with a mean 
measurement error 2.9 km s~^. This is similar to that of the CG clouds (Nielsen et al. 

1998) of 5.1 - 7.4 km s~^. This strongly supports the association of the PMS stars with the 
CG 30/31/38 clouds. The 2.7 km s~^Ylsr of XRS 10 is not significantly lower than the rest 
of the PMS stars. The radial velocities of the stars along the lower locus are inconsistent 
with Vlsr of CGs, unless some of these are single-lined spectroscopic binaries. 

3.6. Spectroscopic Bineiries 

Binary stars yield fundamental astrophysical parameters: stellar masses and occasion- 
ally, radii. It is very difficult to estimate masses and ages of PMS stars by only fitting 
models, as it is the case in our study. Dynamical determination of masses of binary PMS 
stars would enable us to better estimate the ages. 

Spectroscopic binaries have short (few days) periods. This means that the separation 
of two stars are very small. Mayor & Mermilliod (1984) suggested that the dynamical 
evolution of the short period low mass binary systems depends on their PMS history. These 
spectroscopic binaries are, therefore, a very important resource to help us better understand 
the star formation history of multiple systems. 

We identified two double-lined spectroscopic binary systems (SB2): KWW 1637 and 
KWW 1125 (in Figures 5b, 5d). The separation between the Li I 6708 A lines is ~0.7 - 0.8 
A (~33 km s-i) for XRS 6, and -4.5 A (~200 km s"^) for KWW 1125. XRS 6 is a strong 
X-ray source which also shows significant X-ray variability. KWW 1125 has broad Ha and 
Li absorption lines with a large separation between the lines of the two components, which 
means that this system may have a very short, ~few days, orbital period. 

HBC 554, also known as PHa 14 (Pettersson 1987), is a visual binary star (Reipurth 
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& Zinnecker 1993). Therefore, of the 17 stars hsted in the Table 2, there are at least three 
binary systems. Unfortunately, we have only single epoch observations on these systems; we 
can not determine any orbital parameters with only the current data. 

3.7. Color-Color Diagrams and the Circumstellar Disk Fraction 

In the color-color diagrams (Figures 7, 8) we plot the young stars using filled circles (K 
- M stars from the upper locus) and triangles (lower locus) . In Figure 7 the X-ray sources 
are plotted as crosses and asterisks for comparison. 

Among the stars in the upper locus, only KWW 873 shows a convincing near-IR excess in 
both Figures 7 and 8, with H - K ~ 0.2 mag. KWW 1892 and KWW 1637 may have marginal 
excesses. KWW 1892 is clearly a CTTS, based on its Ha equivalent width. However, a 
small excess (H - K <0.1) of KWW 1892 is within the uncertainties of spectral type and 
Ay measurements. KWW 1637, a SB2 binary system, shows variability in X-ray, and also 
between the near-lR photometry from this study and 2MASS photometry, therefore it is 
unclear if this system has ncar-IR excess (H - K ~0.1) due to a circumstellar disk. Four 
upper locus stars lie along the CTTS star locus (Meyer et al. 1997) in Figure 8, but they 
have reddened photospheres with Ay values between 0.5 - 1.3, and do not have near-IR 
excess. The stars in the lower locus (triangles) have normal photospheric colors within the 
scatter. 

A near-IR excess may indicate the presence of an inner circumstellar disk around a 
young star, e.g., (Hillenbrand 2003, and references therein). In our color-color diagrams, 
one PMS star among the eleven on the upper locus shows a convinving near-IR excess. 
The circumstellar disk fraction for this group of stars is ~9%. The disk fraction derived 
using only JHK colors will under-cstimate the true disk fraction, especially among the cooler 
stars, therefore we consider the ~ 9% as a lower limit. This low circumstellar disk fraction is 
consistent with those seen in star forming regions where O stars are present, such as Upper- 
Sco (Walter et al. 1994; Hillenbrand 2003) or NGC 2362, which shows about a 10% disk 
fraction at the age of ~ 5 ±1 Myr. However it is lower than some other known star forming 
regions, such as NGC 2264 (~55% at ~ 3 Myr, Haisch et al. 2001), and the TW Hydrae 
association (~ 20% at ~ 10 Myr). 

The survival time of inner circumstellar disks in T associations is suggested to be about 
3-10 Myr, however exceptions exist in some cases which show very low circumstellar disk 
fraction in young (<3 Myr) associations, e.g., (Walter et al. 1988; Haisch 2001; Hillenbrand 
2003, and references therein). When exposed to strong UV radiation from nearby O stars. 
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surrounding reservoir gas that supplies gas to the circumstellar disks evaporate quickly. By 
losing reservoir gas due to the photoevaporation, those circumstellar disks exposed to UV 
radiation would likely have shorter lifetime than the systems without a source of UV radiation 
nearby. 

4. Discussion 

We have identified 16 young low mass stars along the hne of sight to the CG 30/31/38 
complex. These stars trace out two distinct loci in the CMD. 

The eleven stars that lie along the upper locus have K - M spectral types. Their Li 
abundances are consistent with an age < 5 Myr. Their location in the CMD is consistent 
with a 5 Myr isochrone at the ~200 pc distance of the GC complex. The 4.2 to 8.9 km s~^ 

(±2.9 km s~^) stellar radial velocities agree well with the 5 - 7.5 km s~^ radial velocities of the 
CG 30/31/38 complex (Sridharan 1992; Nielsen et al. 1998). This supports the interpretation 
that these eleven PMS stars are physically associated with the CG complex. 

The stars on the lower locus must be in the background, whether or not they form a 
coherent group. However, they show neither coherent radial velocities (see Table 2) nor a 
spatial concentration. We suggest that these stars are unlikely to be associated with CGs 
or with the 11 PMS stars along the upper locus. These stars all have spectral types earlier 
than K and all lack Ha emission. Five of the F - G stars exhibit strong Li I absorption 
hues. Of the three X-ray sources on the lower locus (KWW 1055, XRS 11, and XRS 12) 
only KWW 1055 shows Li I A 6708 A hne. 

KWW 314, which lies on the lower locus, is an early- type star. We obtained a low 
dispersion (4.3/AA/ resolution) blue spectrum using the SMARTS 1.5m RC spectrograph in 
February 2004. The spectral type is A3; the extinction Eb-v is 0.9 ±0.2 mag. There is no 
emission above the continuum between 3500 and 5100A. The (B — V) — 1.11 is consistent 
with a reddened early A star reddened by Ay ~ 3 mag. This suggests an intrinsic V 12. 
The distance modulus then places this star at ~1.2 kpc. The extinctions to the other stars 
on the lower locus are much smaller, and inconsistent with that of KWW 314. We conclude 
that the lower locus is not a physical association, but a spurious alignment of 6 young, 
possibly PMS, background stars and a perhaps unrelated A3e star. 

The \og{fx/ fboi) values of the PMS stars imply that all the K and M stars detected 
in X-rays are active PMS stars. Since the X-ray detections are only slightly above the 
sensitivity limit, we are likely seeing only the peak of the luminosity distribution if the 
^og{fx/ fboi) relation of this region is similar to that of other star forming regions, e.g., 
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Flaccomio et al. (2003a,b). If we assume that the log{fx/fboi) of this star forming region 
follows the distribution of Orion Nebula Cluster (ONC) (Flaccomio et al. 2003a,b), in the 
mass range of O.2M0 - I.OMq, there should be roughly six times as many PMS stars, all with 
^og{fx/ fboi) < -3.0, as we have detected here. This suggests that deeper X-ray exposures 
may find about 35 more PMS stars with log{fx/ fboi) < -3.0 in this region. 

The V.V — I CMD (Figure 6) shows a cut-off of the upper locus at V ^ 18, below 
which no stars are found. The optical photometry is complete to V ~ 20, therefore the lack 
of stars below V ^ 18 (M ~ 0.2 Mq) is significant. This may be due to selective reddening 
of lower mass stars, or it may be interpreted as an evidence of a truncated mass function 
below a certain mass. Deeper and more complete X-ray and IR observations will help us to 
obtain more complete census of PMS stars to answer whether or not the mass function of 
this star forming region is indeed truncated, and differ from other star forming regions. 

In the color-color diagrams (Figures 7, 8), about 9% of the PMS stars show the near- 
IR excesses indicative of an inner circumstellar disk. This low circumstellar disk fraction, 
also seen in other OB associations, may be due to the influence of ionizing sources, such as 
C Pup, the progenitor of Vela SNR, and other O stars. The parent gas cloud, which feeds 
the accretion disks, may be destroyed by the photoevaporation process from nearby O stars 
at the early stages of PMS evolution. This photoevaporation could result in a short disk 
survival time as well as low final masses of those stars. The lifetime of truncated disks with a 
mass of 1 Mq is about 1 Myr if the mass accretion rate onto the star is about 10~^ Mq yr~^. 
Therefore, about 1 Myr after the disk ceases to accrete from the surrounding cloud the 
circumstellar disks should dissipate, which may explain the low circumstellar disk fraction 
seen here compared to that of quiet star forming regions like Taurus. 

Three PMS stars show evidence of accretion in their Ha emission line profiles. Even 
though their near-lR excesses are not significant, they may still possess circumstellar disks. 
We suggest that the circumstellar disk fraction and disk lifetime depends not only on age, 
but also the environment. Mid-IR photometry will help us to derive a more reliable estimate 
of the circumstellar disk fraction of these systems. 

We identified two double-lined spectroscopic PMS binary systems (KWW 1637 and 
KWW 1125). KWW 1953 was observed twice on different nights: the Li I A6708 line profiles 
are very wide (about 2A) and are also varying. This may due to the fast rotation of the star, 
or could also be a potential double-lined binary that were not completely resolved. 
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4.1. Star formation in photoevaporating globules &; Star formation history of 

CG 30/31/38 complex 

For globule-size triggered star formation, the triggering mechanism can be shock com- 
pression by a supernova explosion (Vanhala at al. 1996), an expanding H II region (Fukuda 
et al. 2000), radiation driven implosion (Lefloch & Lazaref 1994, 1995; Lefloch et al. 1997), 
and the stellar wind from massive stars (Boss 1995, and references therein). Ionizing pho- 
tons from OB stars photoevaporate the molecular clouds, destroying giant molecular clouds. 
However before the total destruction of whole cloud occurs, pre-existing dense parts or cores 
of the parent molecular cloud can survive and form one star or a group of stars. 

The CG formation scenario was discussed by Bertoldi & McKee (1990). Lefloch & 
Lazaref (1994) imply that these small globules in H II regions are evaporating and being 

compressed due to the external ionizing sources. Theoretical studies involving "rocket ef- 
fects" and the "radiation driven implosion" models (Oort & Spitzcr 1955; Rcipurth 1983; 
Lefloch & Lazaref 1994) explain how a shock can drive into a spherical clump, compress it, 
and produce the comctary shape of clouds in pressure equilibrium with surrounding pho- 
toevaporating gas. The globule then can be accelerated away from the UV source as it 
photoevaporates. 

Whether or not small scale triggered star formation can happen in globules depends 
on a number of physical conditions, such as the time scale needed to form a star, the rate 
of photoevaporation, the shock velocity, magnetic field strength and orientation, etc. The 
theoretical studies mentioned above show that star formation can be triggered in globules 
on a time scale short compared to the time needed to destroy or evaporate globules via UV 
radiation and/or shocks. 

Here we suggest a simple star formation history of the CG 30/31/38 complex in the 
broader context of the star formation history of the Gum Nebula. 

About 6-15 Myr ago, young clusters (now Vela 0B2 and Tr 10) formed in a giant 
molecular cloud; OB stars in these clusters ionized their neighboring clouds and the H II 
region expanded (the Gum Nebula). Their wind-blown shells may have been comparable in 
size to those commonly seen in the giant H II regions of other galaxies, such as those seen 
in the HST image of NGC 4214. 

Approximately 5 Myr ago, the clouds photoevaporated in the UV radiation field of 
the luminous O stars, while low mass stars formed in dense cores of the pre-CG 30/31/38 
complex. Under the strong UV radiation from the O stars, the surrounding gas reservoir 
of these young stellar objects began to dissipate. This resulted in accretion disk lifetime 
shorter than those of stars forming far from sources of ionizing radiation. 
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Some 1 million years ago the runaway early O star ( Pup moved closer to this CG 
complex, dissipating the clouds surrounding of these young stars. The kink in the direction 
of the CG 31 tails is consistent with direction of motion of ( Pup. By now, the UV radiation 
has evaporated most of the clouds, though a few dense cores remain. The most recent 
episode of triggered star formation has occurred in the dense core of CG 30, as evidenced by 
the Herbig-Haro object HH 120. Note that the average age of this kind of HH objects are 
suggested to be <10^ years, e.g., (Reipurth et al. 2000). 

5. Summciry 

1. Using multi- wavelength photometry and spectroscopy, we identify 14 new PMS stars, 
adding to the three previously known Ha -emitting stars in the CG 30/31/38 region. The 
spectroscopic sample is complete to V~16.5; the optical limiting magnitude is about 20. 

2. The stars in the CG 30/31/38 complex lie along two distinct loci in the CMD. 

3. The upper locus has an age of < 5 Myr at d = 200 pc. The spectral types of the 
stars on the upper locus are K6 - M4. The mean Li abundance for M stars, logN(Li) = 2.68 
±0.28 for NLTE and 3.26 ±0.15 for LTE), implies an age of 2 - 5 Myr and supports the 
youth of the PMS stars. The radial velocities of these stars are consistent with the published 
velocities of CG 30/31/38 cloud complex. 

4. The lower locus has an age of < 100 Myr at d~ 2 kpc. The stars along the lower 
locus have spectral types F - G, with one Ae star. With an exception of KWW 1055, these 
stars have radial velocities that are not consistent with those of the CG 30/31/38 complex. 

5. The ~9% circumstellar disk fraction derived using optical and near-IR excesses is 
consistent with those seen in other star forming regions with O stars, such as the Upper-Sco, 
and is lower than those typical of quiet regions, such as Taurus. This suggests that the loss 
of the gas reservoir for accretion disks at this early stage of star formation is due to the 
radiation from the OB stars. 

6. We find two double-lined spectroscopic binaries in the sample. Including the known 
visual binary HBC 554, the observed multiplicity fraction among the 17 PMS stars is 18%, 
which sets the lower limit to the true multiplicity fraction in this region. 

8. We conclude that there have been at least two episodes of star formation in the 
CG 30/31/38 region: (1) On-going star formation triggered by UV radiation from OB stars 
in the head of CG 30 cloud, as exemplified by HH 120. (2) < 5 Myr old low mass (0.2 - 1 
M0) PMS stars that outline the CG 31 complex, whose formation may have been triggered 
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by pre-existing O stars like the progenitor of Vela SNR and ( Pup, as well as other OB stars 
in Vela 0B2 and Tr 10. 

We thank Dr. WiUiam Sherry at NSO/NOAO and Dr. Rob Jeffries at Keele University 
for helpful discussions. We also thank CTI0/4m crew who helped us to use HYDRA multi- 
fiber spectrograph. This research was supported in part by NASA grant NAGS- 1594 to 
SUNY Stony Brook. 

A. Notes on individual PMS stars 

In this appendix we discuss some of the spectroscopically confirmed PMS stars. Coor- 
dinates, optical, and near-IR photometry are in Tables 3 and 4 in the discussion section. 

KWW 464 

This star 17" away from XRS 1, a 4 cr deviation. The M3 spectral type is based on the TiO 
band strengths. VFx(Li I) ~0.66A indicates that this is a young star that is unlikely to be 
older than 5 - 10 Myr. The equivalent width of Ha is ~-4A, therefore it is a NTTS/WTTS 
(see Table 2). 

KWW 1892 = XRS 2 = HBC 553 

Figure 5a shows two spectra of KWW 1892 (XRS 2; Ml - M2), obtained on consecutive 
nights. This object, also known as HBC 553 (Pettersson 1987; Herbig & Bell 1988) has the 
strongest Ha emission line {Who ~-30A) in our sample. The Ha emission line profile shows 
changes on the second night, while the absorption line profiles do not. The FWHM of Ha 
line varied from 4 - 5A (night 1) to 3A (night 2). Based on Who, this is the only classical 
T Tauri star in our sample. 

KWW 1637 = XRS 6 

Li I A 6708A, Ca I A 6717A, and other Fe I lines all are doubled, strongly indicating that 
it is a double-lined spectroscopic PMS binary system which has 0.8 - lA separation (40 - 
45 km s~^) at this epoch. The X-ray emission was strong and variable. 

KWW 873 = XRS 7 

Although Whu is less than lOA, the line profile, with a blue shifted absorption reversal, 
implies a strong wind in this system. 

KWW 975 = HBC 554 

The spectrum of KWW 975 (HBC 554) in Figure 5b shows asymmetric Ha and Li 1 line 
profiles. The Ha line profile has dips (stronger on the redshifted side), which could indicate 
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accretion or infalling gas in the system. Background sky emission subtraction affects only 
tlie narrower component, and not tlie broader line profile. Herbig & Bell (1988) identified 
HBC 554 as a visual binary system, and classified the spectral type of this star as M 1.5, 
which is consistent with our spectral type estimate (M 2). 

KWW 1953 

The Li I 6708A line and Fe I lines show asymmetric and broad line profiles. KWW 1953 
may be a rapidly rotating star, or a unresolved double line binary system. 

KWW 1055 = XRS 10 

The Ha line is in absorption (Figure 5d). The Li line equivalent width is consistent with 
no Li depletion. Although appear on the lower locus, and hence at a greater distance, the 
radial velocity is marginally consistent with that of the CG 30/31/38 complex. 

XRS 9 

XRS 9 is a fast rotator. It is not a Ha emission source, but shows a marginal Li I line 
(Figure 5d). All the absorption lines are broadened, most likely due to rotation. All the 
lines, including the metal lines Fe I and Ca 1, show fiat bottoms which may suggest a spot 
on its surface. The radial velocity is inconsistent with those of the of CGs and associated 
PMS stars. It is likely a young star in the background. 

KWW 314 

The double-peaked Ha emission arising from a broad Ha absorption line is typical of a Ae/Be 
star (and not necessarily of a young object). A low dispersion spectrum confirms that this 
is an A3 (±2) star with Eb-v = 0.9±0.2 mag. A^, is '^3 mag, (see tables 3 and 4), and 
thus V ~12mag. The distance modulus is then ~ 10.5 magnitude (~1.2 kpc). The radial 
velocity {Vlsr — -15.11 km s~^) is inconsistent with that of the CGs and PMS stars. This 
is therefore likely a background star unrelated to the CG 30/31/38 complex. 

KWW 1125 

This is a young double-lined spectroscopic binary system with undepleted Li abundances. 
The spectrum (Figure 5d) shows broad and double peaked Ha and Li I absorption lines. Fe I 
and Ca I fines also seem to have double fines. The two components are ~4A (180 km s~^) 

apart at the epoch of observation; this is likely to be a short period system with a period 
of a few days. The redshifted component has ~50% of the line strength of the blueshifted 
component in Ha , but the Li I strengths of both components are similar. At this epoch the 
radial velocities are -65.6 km s~^ for the blue shifted component and 152.6 km s^^ for the 
redshifted component. Since only one epoch was observed, we have no information on the 
systemic velocity or the mass ratio. 
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Table 1. ROSAT/HRI sources in CG 30/31/38 complex 
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Note. — Values are taken from the Standard Analysis Software System (BASS) output. Note that the 
uncertainty on the count rate docs not directly translate to a Signal to Noise (S/N) ratio: the former 
corresponds to a larger radius containing 90% of the encircled energy. The significance of the count rate 
of the weaJf sources is decreased because background noise dominates. 

"■X-ray - optical position offset in RA 

''X-ray - optical position offset in DEC 

■^Detect cell size in arcsec 

"^Net count rate within the 90% encircled energy radius 

''Net source counts (above background) within the 90% encircled energy radius 
'^signal-to-noise ratio within the detect cell 
^Not detected in the 24 arcsec detection cell 



Table 2. Young stars identified from this spectroscopic study 
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Note. — The spectral types and equivalent widths of Ha A6563A and Li I A6708A for PMS stars are presented. Note that XRS 9, KWW 1055, KWW 1125, KWW 1333, 
and KWW 1806 are more likely 50 - 100 Myr old field stars with strong magnetic activity. 

''The errors were determined by measuring W;^(Li I) multiple (>10) times. 

t'ithium abundance for both LTE and NLTE model (NLTE from Pavlenko & Magazzu 1996; LTE from Duncan 1991) 

'^The uncertainty in Ay value is mainly introduced by the uncertainty of spectral type determination (~0.5 - 1 subtype for K and M stars, >1 for G type and earlier). 

dPettersson (1987), Hcrbig & Bell (1988) 

^Both components show the same Lithium strength. The equivalent width here is for the blueshifted component. 
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Table 3: Optical photometry 
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Table 4: Near IR photometry 



ID K aK J - K a{J - K) H - K a{H - K) J° crj W uH K°- aK 
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a: JHK data from the 2MASS All-Sky Release Point Source Catalog. 
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Fig. 1. — The DSS image of the CG30/31/38 complex, showing the HRI sources and the 
PMS stars. North is up; east is to the left. CG 30 has an outflow source HH 120; CG 31 is 

divided into 5 pieces, and CG38 is the smallest globule, to the southeast of CG 30. The ten 
ROSAT HRI sources (boxes) appear to surround the heads of the CGs on the side toward 
the central ionizing sources, while the CG tails point away from the ionizing sources. The 
small circles are centered on the spectroscopically confirmed PMS stars that arc associated 
with CGs, while larger circles are stars that are less likely to be associated with CGs. 
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Fig. 2.— The V,V - I color magnitude diagram (CMD) of CG 30/31/38 complex. Stars 
which are likely counterparts of X-ray sources are indicated with crosses and asterisks. Di- 
amonds mark Ha emission sources (Pettersson 1987; Herbig & Bell 1988; Schwartz et al. 
1990). The Ay=l mag reddening vector, plotted at top center, runs almost parallel to the 
isochrones. Isochrones are from Baraffe et al. (1998) using color tables by Kenyon & Har- 
mann (1995); Luhman (1999); Leggett et al. (2001) for 2 Myr, 5 Myr, and 2 Gyr at an 
assumed distance of 200 pc. We use the evolutionary models with [M/H]=0 and y=0.275. 
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Fig. 3. — The log fx, log/(,o; plot for ROSAT/HRI sources. The dotted hne is drawn at 
^og{fx/fboi) = -3. 



-36- 



o 
J2 



-1.0 
-1.5 

-2.0 

h2< -2.5 

-3.0 

-3.5 
-4.0 



o 



3.8 



3.7 3.6 

log(Teff) 



3.5 



Fig. 4. — ^og{fx/ fboi) plotted as a function of Teff- The dashed hne shows the approximate 
3(7 detection hmit. Half of the X-ray sources are on or slightly over the sensitivity limit, 
which suggests that there are more PMS stars awaiting to be discovered. 
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Fig. 5a. — HYDRA/echelle spectra of confirmed PMS stars associated with cometary glob- 
ules. Continuua are flattened and normalized to unity. 
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HYDRA/echelle spectra of confirmed PMS stars associated with cometary glob- 
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HYDRA/echelle spectra of confirmed PMS stars associated with cometary glob- 
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Fig. 5d. — HYDRA/echelle spectra of young stars that are unlikely to be associated with 
the CG 30/31/38 complex. 
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Fig. 5e. — HYDRA/echelle spectra of young stars that are unlikely to be associated with 
the CG 30/31/38 complex. 
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Fig. 6. — Spectroscopically confirmed PMS stars in the V,V — I CMD. This is basically 
the same plot as Figure 2. Filled circles and open triangles are spectroscopically confirmed 
PMS stars including KWW 1043 (HBC 555) along the upper PMS locus, and the A3e star 
(KWW 314) along the lower locus. The stars in the upper PMS locus (filled circles) are 
mainly K and M stars, and appear under the head of CG 31. The stars in the lower locus 
(open triangles) are F - G type stars, and appear primarily around the CG 30. The asterisks 
are optically-unobserved X-ray sources. The isochrones are the same as in Figure 2. 
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Fig. 7. — The I — K,V — I color-color diagram. Filled circles mark spectroscopically- 
confirmed PMS stars that are kinematically associated with the CGs. The open triangles 
mark stars that lie along the lower locus, including the A3e star KWW 314, and open 
diamond symbols are Ha sources. The isochrones (Baraffc ct al. 1998, Luhman 1999) arc 
plotted for 2 Myr and 2 Gyr and distance of 200 pc. The reddening vector is shown upper 
right hand side of the panel for Ay = 1, and the dash-dot lines indicate reddening directions. 
The label "IR excess" indicates the region of the plot occupied by IR excess objects. One 
source KWW 873 shows I - K excess. 
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Fig. 8. — The J — H, H — K color-color diagram of the field. The symbols used are the 
same as in Figure 7. The dash-dot fines are the CTTs locus from Meyer et al. 1997, and the 
solid lines indicate approximate locations of unreddened dwarf stars using color tables by 
(Kenyon & Harmann 1995) modified using (Luhman 1999; Leggett et al. 2001) as in Table 3 
in (Sherry et al. 2004). 



